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Abstract 

For any MONDian extended theory of gravity where the rotation curves of spiral galaxies are explained 
through a change in physics rather than the hypothesis of dark matter, a generic dynamical behavior 
is expected for pressure supported systems: an outer flattening of the velocity dispersion profile 
occurring at a characteristic radius, where both the amplitude of this flat velocity dispersion and 
the radius at which it appears are predicted to show distinct scalings with the total mass of the 
system. By carefully analyzing the dynamics of globular clusters and elliptical galaxies, we are able to 
significantly extend the astronomical diversity of objects in which MONDian gravity has been tested, 
from spiral galaxies, to the much larger mass range covered by pressure supported systems. We show 
that a universal projected velocity dispersion profile accurately describes various classes of pressure 
supported systems, and further, that the expectations of extended gravity are met, across seven orders 
of magnitude in mass. These observed scalings are not expected under dark matter cosmology, and 
would require particular explanations tuned at the scales of each distinct astrophysical system. 

Keywords: gravitation — stars:kinematics and dynamics — galaxies:star clusters:general — galaxies: 
fundamental parameters — galaxies: kinematics and dynamics 


1. INTRODUCTION 

Within the context of modified gravity theories con¬ 
structed to yield galactic dynamics in the absence of 
dark matter, one of the most convincing supporting ar¬ 
guments comes from the success of MOND in repro¬ 
ducing rotation curves of spiral galaxies (Milgrom 1983; 
Swaters et al. 2010; McGaugh 2012). Both the flatness 
of the observed spiral rotation curves and the scaling 
of the asymptotic rotation with the fourth root of the 
total baryonic content, i. e. the Tully-Fisher relation, 
can be accurately reproduced across spiral galaxy sizes 
and types, as well as the overall shapes and correlations 
between observed features in the baryonic distribution 
and small kinematic irregularities (Sanders & McGaugh 
2002; Famaey & McGaugh 2012; Rodrigues 2014; Lelli 
et al. 2016). However, such successful tests are re¬ 
stricted to the range of mass scales and morphology of 
spiral galaxies. 

To test the validity (or lack thereof) of MONDian 
gravity schemes, it is desirable to extend the diversity 
of astrophysical systems over which these ideas can be 
probed. To do this, we give the expected kinematic pro¬ 
files for two different types of pressure supported sys¬ 
tems, astrophysical objects with baryonic masses span¬ 


ning from IO^Mq for globular clusters, to IO^^Mq for 
elliptical galaxies. 

By the term MONDian gravity we shall refer to any 
extended gravity theory where at a > oq scales Newto¬ 
nian gravity is recovered, while the a < oq regime re¬ 
produces MOND dynamics, at the v << c limit, where 
qq is Milgrom’s acceleration of 1.2 x 10 “^°to/s^ (Beken- 
stein 2004; Gapozziello & De Laurentis 2011; Mendoza 
et al. 2013; Verlinde 2016). Generically for such the¬ 
ories, beyond a radius given by Rm = 
centrifugal equilibrium velocities will become flat at a 
level of U = (GMoo)^'^^ (Milgrom 1983). Similarly, for 
pressure supported systems, beyond around Rm, veloc¬ 
ity dispersion profiles will stop falling along Newtonian 
expectations and flatten out at a level of CToo = U/3^/^, 
with M the total baryonic mass of a rapidly converg¬ 
ing astrophysical system (Milgrom 1984; Hernandez & 
Jimenez 2012). 

Recently, McGaugh et al. (2016) and Lelli et al. 
(2016) used the SPARG rotation curve and photome¬ 
try databese, to show that the baryonic mass distribu¬ 
tion alone directly determines the resulting acceleration, 
along Newtonian expectations in the high acceleration 
regime, and along MONDian ones when a < qq. A 
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similar study reaching conclusions also in support of 
MONDian dynamics, but usiirg structural, dyrramical 
and mass determinations for a range of elliptical galax¬ 
ies is that of Dabringhausen et al. (2016). Baryonic 
mass determination in the above studies however, are 
subject to numerous assumptions, as well as uncertairr- 
ties and systematics, which vary for differirrg classes of 
astrophysical systems. Hard to determine gas fractions, 
star formation histories, the presence of multiple and 
complex stellar populations and unknown stellar mass 
functioirs imply mass to light ratios and their radial vari- 
atioirs are subject to large empirical errors. 

Here we take a complementary approach, using only 
kinematical data and avoiding any dynamical physical 
assumptions or mass estimates in deriving the param¬ 
eters we study, concentrating merely on the shape of 
observed projected velocity dispersion profiles. The ex¬ 
pected theoretical scalings of both Rm and CToo with M 
imply the mass can be eliminated to obtain Rm oc cr^, 
a relation between the parameters of the directly ob¬ 
servable velocity dispersion profiles which can now be 
directly tested. 

By analysing velocity dispersion profiles from Galac¬ 
tic globular clusters recently observed by Scarpa et al. 
(2007a; 2007b), Scarpa & Falomo (2010) and Scarpa et 
al. (2011), henceforth the Scarpa et al. group, and Lane 
et al. (2009), Lane et al. (2010a; 2007b) and Lane et al. 
(2011), henceforth the Lane et al. group, and elliptical 
galaxies having high quality 2D kinematics by Sanchez 
et al. (2016b), using data from the CALIFA project, 
(Sanchez et al. 2012; Walcher et al. 2014; Sanchez et 
al. 2016a), we show that a universal projected velocity 
dispersion profile cr{R) = + CToo adequately 

serves to model the two distinct classes of astrophysical 
systems. Further, a scaling of Rm oc cr^ ensues, with 
n = 1.89 ± 0.32, compatible with the generic expecta¬ 
tions of MONDian gravity models of n = 2. 

In section (2) we present the development of simple 
first order predictions for the large scale velocity dis¬ 
persion expected under MONDian gravity, and the cor¬ 
responding scalings between the parameters describing 
such velocity dispersion profiles for pressure supported 
systems. Section (3) includes a description of the data 
used for sets of velocity dispersion profiles for globular 
clusters and elliptical galaxies, together with the fitting 
procedure and parameters obtained for the universal ve¬ 
locity dispersion profile proposed. In section (4) we show 
how the two sets of systems treated, spanning seven or¬ 
ders of magnitude in mass, follow the proposed univer¬ 
sal velocity dispersion profile, and show also how the 
flattening radii follow an overall scaling consistent with 
being proportional to the square of the asymptotic veloc¬ 
ity dispersion, as predicted generically under MONDian 
gravity. Section (5) summarizes our conclusions. 


2. MONDIAN THEORETICAL EXPECTATIONS 

Thinking of modified gravity in terms of a Newtonian 
force law implies a change of regime from the Newtonian 
expression of Fjv = GM/r'^ to Fm = {GMao)^/'^/r, oc¬ 
curring at scales of Rm = (GM/oq)^/^, if one wants 
to model observed galactic dynamics in the absence of 
any dark matter (Milgrom 1983). The details of the 
transition are not clear, beyond the requirements from 
consistency with solar system dynamics (Mendoza et al. 
2011), Milky Way rotation curve comparisons (Eamaey 
& Binney 2005) and Galactic globular cluster modelling 
(Hernandez et al. 2012), for the transition to be fairly 
abrupt. Under such schemes, centrifugal equilibrium ve¬ 
locities will tend to the Tully-Eisher value of: 

U=(GMao)'/^ (1) 

for test particles orbiting a total baryonic mass M. Un¬ 
der the MONDian regime, the relation between centrifu¬ 
gal equilibrium velocities and isotropic velocity disper¬ 
sion will be only slightly modified with respect to the 
Newtonian case, to yield: 

aoo = V/Vi, (2) 

e.g. Milgrom (1984), Hernandez & Jimenez (2012). 

Thus, we can write the isotropic Maxwellian velocity 
dispersion of a pressure supported system in the MON¬ 
Dian regime as: 

<jI = ^ {GMa^f'^ . (3) 

In order to extend the range over which MONDian dy¬ 
namics can be tested to the much larger ones over which 
pressure supported systems appear, would require ob¬ 
servations of CToo and independent inferences of the total 
baryonic mass of various systems, to directly test the 
validity of Equation (3). Any such undertaking will be 
hindered by the large theoretical uncertainties in total 
(and radially changing) mass to light ratios, uncertain 
gas fractions, and details of the relevant IMEs and star 
formation histories, which to make matters worse, vary 
significantly from the “simple” cases of globular clusters 
to the more complex histories of elliptical galaxies. One 
can circumvent such uncertainties and obtain a simpler 
empirical test by eliminating M from the above equation 
in favor of Rm to yield: 

Rm = (4) 

Oo 

which in astrophysical units reads: 
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The above expression now relates only features of the 
kinematic velocity dispersion profile of an astronomical 
system, as aoo will be the large scale asymptotic pro¬ 
jected velocity dispersion profile, and Rm will be the 
radial scale at which the velocity dispersion profile be¬ 
comes flat. The above will only make any sense if ob¬ 
served velocity dispersion profiles do indeed show such 
morphology, for a decline in an inner Newtonian region, 
and a transition to a flat velocity dispersion regime on 
crossing the Rm threshold. 

We are inspired to attempt such a comparison across 
the range of scales covered by pressure supported sys¬ 
tems by recent results coming from the observed ve¬ 
locity dispersion profiles of Galactic globular clusters 
(Scarpa et al. 2010; Scarpa & Falomo 2011) showing 
precisely such a morphology, and recently shown to com¬ 
ply with Tully-Fisher MONDian expectations (Hernan¬ 
dez et al. 2013a). Further, Desmond (2017) recently 
demonstrated through an extensive sample of rotation 
curve observations and baryonic mass inferences across 
sizes and galactic types, that the inferred mass discrep¬ 
ancy acceleration relation is consistent with MONDian 
expectations. 

In the following section we show a number of veloc¬ 
ity dispersion profiles having precisely the morphology 
described above, for Galactic globular clusters and el¬ 
liptical galaxies, which in Section 4 we compare to the 
expectations of Equation 5. 

Notice also that multiplying Equations (3) and (4), 
and substituting M on the right hand side for CToo, the 
resulting combination of parameters Rm<^^ = GM/2> is 
highly reminiscent of the corresponding Newtonian ex¬ 
pression, changing Rm for the half mass radius, and Uoo 
for the central velocity dispersion of a purely Newtonian 
system (Gappellari et al. 2006; Wolf et al. 2010). 

In the above we have assumed that the total baryonic 
mass has essentially converged, and on reaching the flat 
dispersion velocity region the potential of the system 
in question can be treated as a point mass, in consis¬ 
tency with the validity of Newton’s theorems for spher¬ 
ically symmetric mass distributions also under MON¬ 
Dian gravity (e.g. Mendoza et al. 2011), and that the 
dynamical tracers from which Uao is calculated behave as 
test particles. Although this will never be strictly true, 
it will hold to a good approximation, given the very cen¬ 
trally concentrated mass profiles of e.g. Galactic globu¬ 
lar clusters and the inferred density profiles of close to 
the de Vaucouleurs surface brightness of elliptical galax¬ 
ies. Indeed, pressure supported systems have been accu¬ 
rately modeled recently in the MONDian regime repro¬ 
ducing simultaneously the observed surface brightness 
profiles and the measured velocity dispersion ones, e.g. 
Galactic globular clusters by Gentile et al. (2010), the 
giant elliptical galaxy NGC4649 by Jimenez et al. (2013) 


and tenuous galactic stellar halos by Hernandez et al. 
(2013b), all in absence of any otherwise dominant dark 
matter component. More recently, Ghae & Gong (2015) 
studied the predicted distribution of velocity dispersion 
profiles (VDP) for elliptical galaxies and found that their 
MONDian models could reproduce the observed distri¬ 
bution of VDP slopes requiring only the observed stellar 
mass, without any dark matter, and Tian & Ko (2016) 
concluded that the dynamics of seven elliptical galax¬ 
ies, as traced by stellar and planetary nebulae, could be 
well explained by MOND. The issue remains controver¬ 
sial, as some studies have found dark matter halos to be 
a better fit than MOND in elliptical galaxies (Richtler 
et al. 2008; Samurovic 2014, 2016). 

3. EMPIRIGAL VELOCITY DISPERSION 
PROFILES 

3.1. Galactic globular clusters 

In data recently available from the Scarpa et al. group 
and the Lane et al. group, not only the central values, 
but the radial projected velocity dispersion profiles for 
several systems have been published, extending out to 
several half-light radii. As first mentioned by the Scarpa 
et al. group, a flattening of the velocity dispersion pro¬ 
files for Galactic globular clusters appears on crossing 
the oq acceleration threshold, in accordance with MON¬ 
Dian gravity expectations. In Hernandez & Jimenez 
(2012) some of us recently modeled the above globu¬ 
lar cluster data to construct self-gravitating models for 
the systems studied, where both the projected surface 
brightness profiles and the observed velocity dispersion 
profiles have been accurately reproduced, using a mod¬ 
ified gravity force law of the type described in section 
(2). Further, in Hernandez et al. (2013a) we showed 
that the asymptotic velocity dispersion values for the 
clusters treated actually match the expected scaling of 
Coo ~ {GMa^y^'^ of MONDian gravity, through com¬ 
paring with stellar synthesis models tuned to the ages 
and metallicities of each of the individual globular clus¬ 
ters studied. In the above works we also showed that 
an empirical projected velocity dispersion profile of the 
form: 

a{R) = -I- (Too (6) 

serves to accurately model the reported velocity disper¬ 
sion profiles of Galactic globular clusters. In Equation 
(6) (Too is the asymptotic value for the system’s velocity 
dispersion, the central value for this quantity is given 
by (t(0) = ao + CToo, and R^ gives a transition radius 
beyond which the asymptote is rapidly approached. 

Here we include again the fits to the data of the Scarpa 
et al. and Lane et al. groups to Equation (6). We 
use a non-linear least squares Levenberg-Marquardt al- 


(R)[km/s] ff(R)[kin/s] ff(R) [km/s] cr(R) [km/s] 
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NCG 6656 


NGC 6752 


NGC 7087 


Figure 1. Projected velocity dispersion profiles for the 12 Galactic globular clusters of our sample, as a function of radial 
distance in the system, with vertical error bars showing the Icr dispersion at each radial bin. The solid curve gives the best fit 
to the universal profile proposed. 
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NGC 0677 





NCG 0938 NGC 0962 NGC 1026 



Figure 2. Projected velocity dispersion profiles for the 13 elliptical of our sample, as a function of radial distance in the system, 
with vertical error bars showing the la dispersion at each radial bin. The solid curve gives the best fit to the universal profile 
proposed. 
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gorithm to estimate the parameters, and obtain self- 
consistent Icr confidence intervals, which are reported 
in Table 1. We have trimmed the total sample used in 
Hernandez et al. (2013a) to 12 globular clusters where 
the relative errors in the fitted parameters are in all 
cases smaller than 50%. This excludes 4 globular clus¬ 
ters with poorly determined velocity dispersion profiles, 
also keeping the globular cluster sample from growing 
much beyond the numbers of available well measured 
systems in the category of low rotation elliptical galax¬ 
ies. This last to ensure that the systematics relevant to 
any particular class of objects do not dominate and bias 
the overall comparisons we attempt. 

Figure 1 provides the observed velocity dispersion pro¬ 
files for the Galactic globular clusters of our sample with 
the best fit to Equation (6) for each system, showing 
clearly an excellent empirical representation of the mea¬ 
sured profiles. A clear decline in the Newtonian inner 
region is evident, followed by a transition to a “Tully- 
Fisher” flat asymptote. 

3.2. Elliptical galaxies 

A sample of low rotation elliptical galaxies was ob¬ 
tained from Sanchez et al. (2016b), in which 200 radial 
2D velocity dispersion profiles were produced from the 
second data release of the CALIFA survey (Sanchez et 
al. 2012; Walcher et al. 2014; Sanchez et al. 2016a) 
using the PIPE 3D tool. To ensure that the selected 
sub-sample was composed of systems with the least dy¬ 
namical support besides velocity dispersion, we selected 
galaxies based on their morphology (E0-E3). Our sam¬ 
ple comprises only extremely slow rotators, with an av¬ 
erage value of a maximum rotation velocity to certral 
sigma per galaxy of Ymaxl^Q = 0.213, very early type 
systems with negligible total gas content. 

We use data directly obtained from Sanchez et al. 
(2016b), consisting of 2D radial stellar velocity disper¬ 
sion observations with their respective uncertainties, as 
well as associated flux for a correct weighted statistical 
development. From the 2D velocity dispersion observa¬ 
tions we extract a projected radial velocity dispersion 
representation. Due to the large number of pointings in 
each galaxy, we averaged the projected velocity disper¬ 
sion observations within 20 radial bins (Figure 2), each 
with a corresponding Icr dispersion which greatly over¬ 
shadows the intrinsic observational error for each data 
point. In spite of the high Icr dispersion towards the 
center of the galaxies, we adopt this method for sta¬ 
tistical consistency with the globular clusters treated. 
Monte Carlo generated errors for each data point are 
only about 15 km/s, while the radial bin average confi¬ 
dence intervals can reach up to 50 km/s. As described 
in Sanchez et al. (2016a,b) velocity dispersion mea¬ 
surements of under 40 km/s from the middle resolu- 



Figure 3. Observed projected velocity dispersion profiles in 
units of the corresponding cjoo values as a function of radius 
normalized to the corresponding values, for two distinct 
systems. Red circles correspond to measurements for the 
Galactic globular cluster NGC 5139, with a central velocity 
dispersion of 15.2fcm/s and Ra = 16.7pc. Green triangles 
correspond to the elliptical galaxy NGC 0677, with a cen¬ 
tral velocity dispersion of 195fcm/s and Ra = 4.3kpc. The 
consistency of the two scaled projected velocity dispersion 
profiles is evident, despite the range of scales covered by the 
comparison. 

tion setup V1200 are not reliable, so we have discarded 
these data points in our estimates. Using again a non¬ 
linear least squares Levenberg-Marquardt algorithm, we 
fit the proposed universal function of Equation (6) to the 
projected velocity dispersion profiles and obtain optimal 
parameters as well as their respective confidence inter¬ 
vals. After discarding various systems with poorly de¬ 
termined parameters, with fractional errors larger than 
50% in the fitted parameters, we compile a final sample 
of 13 elliptical galaxies, as reported in Table 1. 

Figure 2 shows the radial binning of the projected 2D 
radial velocity dispersion of the elliptical galaxies in the 
sample, as well as the best fit to the universal function 
of Equation (6). Note again the suitability of the fits in 
all cases. 

Although the general fit of the proposed universal ve¬ 
locity dispersion profile is clearly very good across the 
range of objects treated (see the full figure set), spe¬ 
cially towards the flattened outer regions, it is also clear 
that in a number of both globular clusters and elliptical 
galaxies, a slight central enhancement in the data with 
respect to the htted function appears. This cuspy ve¬ 
locity dispersion profiles, e.g. the globular cluster NGC 
2419 and elliptical galaxies NGC 0962 and NGC 1060, 
could signal the presence of a central black hole. This 
is a common feature in early type galaxies, and could 
also appear in globular clusters. For example, Baum- 
gardt (2017) recently concluded from dynamical studies 
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Table 1. Parameters for velocity dispersion profiles of the the pressure supported systems. 


Type of system 

Identification 

ao{km/s) 

(Joa(km/s) 

Aa (pc) 

Globular Gluster 

NGC 0104 

4.5 ±0.4 

5.0 ±0.1 

14.8 ±0.9 

Globular Gluster 

NGC 1851 

4.1 ± 1.7 

3.9 ±0.4 

7.9 ±2.3 

Globular Gluster 

NGC 1904 

2.4 ±0.6 

2.0 ±0.4 

10.3 ±3.3 

Globular Gluster 

NGC 2419 

4.7 ±0.8 

1.0 ±0.4 

1.8 ±0.3 

Globular Gluster 

NGC 5139 

8.0 ±0.7 

7.2 ±0.4 

16.7 ± 1.9 

Globular Gluster 

NGC 6171 

1.4 ±0.4 

2.7 ±0.2 

3.8 ± 1.1 

Globular Gluster 

NGC 6218 

3.3 ±0.9 

1.4 ±0.5 

4.7 ± 1.2 

Globular Gluster 

NGC 6341 

3.8 ±0.7 

3.1 ±0.4 

6.9 ± 1.4 

Globular Gluster 

NGC 6656 

3.8 ± 1.3 

3.3 ±0.8 

7.4 ± 1.7 

Globular Gluster 

NGC 6752 

3.5 ±0.8 

2.0 ±0.3 

10.7 ± 1.4 

Globular Gluster 

NGC 7078 

5.1 ± 1.0 

3.0 ±0.3 

8.7 ± 1.7 

Globular Gluster 

NGC 7099 

1.8 ±0.5 

2.2 ±0.3 

6.2 ±2.3 

Elliptical Galaxy 

NGC 0155 

87.6 ±25.5 

69.1 ± 25.4 

(8.4 ±3.6) X 10® 

Elliptical Galaxy 

NGC 0677 

103.6 ± 24.4 

91.4 ± 16.7 

(4.3 ± 1.6) X 10® 

Elliptical Galaxy 

NGC 0731 

129.0 ± 26.9 

41.6 ±6.6 

(1.4 ±0.2) X 10® 

Elliptical Galaxy 

NGC 0938 

124.0 ± 25.8 

59.0 ± 7.0 

(2.9 ±0.5) X 10® 

Elliptical Galaxy 

NGC 0962 

116.3 ± 22.5 

57.3 ±8.5 

(3.7 ±0.8) X 10® 

Elliptical Galaxy 

NGC 1026 

123.6 ±37.7 

98.8 ±8.8 

(1.8 ±0.4) X 10® 

Elliptical Galaxy 

NGC 1060 

165.6 ± 18.8 

79.5 ± 16.7 

(9.7 ±2.0) X 10® 

Elliptical Galaxy 

NGC 4841A 

111.7 ±46.7 

129.9 ±8.9 

(3.2 ± 1.4) X 10® 

Elliptical Galaxy 

NGC 5198 

163.6 ± 27.3 

53.8 ± 14.7 

(1.5 ±0.2) X 10® 

Elliptical Galaxy 

NGC 5216 

105.1 ± 28.3 

50.2 ±9.6 

(1.2 ±0.3) X 10® 

Elliptical Galaxy 

NGC 6515 

85.6 ± 24.3 

88.4 ± 13.5 

(5.8 ±2.2) X 10® 

Elliptical Galaxy 

NGC 7194 

110.4 ±33.6 

132.3 ± 18.7 

(8.0 ±3.2) X 10® 

Elliptical Galaxy 

NGC 7619 

151.3 ± 21.8 

100.3 ± 14.7 

(4.6 ± 0.9) X 10® 


The first two columns give the object class and identifier for the systems studied. The following three entries show the 
parameters of the fits to the observed projected velocity dispersion profiles and their confidence intervals. Data from the 
Scarpa et al. group and the Lane et al. group for Galactic globular clusters, and data from the CALIFA collaboration for 
elliptical galaxies. 


of NGC 5139 that this cluster likely harbors a 40,000 
solar mass black hole in its center. From Figure 1 we 
see that indeed, this particular system shows the most 
prominent central velocity dispersion excess with respect 
to the universal profile fitted of all the globular clusters 
treated. 

Fortunately, this inner velocity dispersion excesses are 
very central features which do not affect the determina¬ 
tions of Ra or tToo j the parameters which can be associ¬ 
ated to MONDian dynamical predictions. 

4. COMPARISONS WITH MONDIAN 
EXPECTATIONS 

Comparing Figures 1 and 2, we see that although the 
physical scale is changing from pc to kpc, the velocity 
profiles of globular clusters and ellipticals can both be 
fitted adequately by the same family of functions. Fur¬ 
ther, while the peak uo is on average different for glob¬ 


ular clusters and ellipticals when scaled with a^c, the 
transition from a declining velocity profile to a flat one 
emerges with remarkable consistency at a scaled radius 
of R/Ra ~ 2. 

From Equation (6) we see that the central velocity 
dispersion values, given by cr(0) = gq + CToo, span two 
orders of magnitude, with the scale radii Rc covering 
three orders of magnitude, and the total baryonic masses 
going from a few IO^Mq scale for the smallest Galactic 
globular clusters to 10 ^^Mq for the elliptical galaxies. 
If we now scale the projected velocity dispersion profiles 
of the various systems presented in the previous section, 
and consider a{R)/aoc as a function of R/Ra for each 
profile, only one degree of freedom remains, given by the 
central velocity dispersion, a(0)/aoo- 

Therefore, if we choose distinct systems having equal 
cr(0)/cr Qo v3,lii6S, they should look idciiticul in n ^oo 

vs. R/Ra plot. Figure 3 shows the scaled velocity dis- 
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Figure 4. Observed projected velocity dispersion profiles in 
units of the corresponding aoo values as a function of radius 
normalized to the corresponding values, for two distinct 
systems. Red circles correspond to measurements for the 
Galactic globular cluster NGC 7078, with a central veloc¬ 
ity dispersion of 8.1km/s and Ra = 8.7pc. Green triangles 
show the elliptical galaxy NGG 6515, with a central velocity 
dispersion of nikmjs and Ra = 5.8kpc. The two scaled 
profiles can be modeled through the same functional form 
of the proposed universal proHle, the Figure also shows the 
tendency for the scaled profiles of elliptical galaxies to be 
flatter and for globular clusters to be steeper. 

persion profiles of two different systems chosen to have 
very similar a{0)/aoo values. Circles correspond to the 
Galactic globular cluster NGC 5139, while triangles rep¬ 
resent data for the elliptical galaxy NGC 0677. The solid 
line shows a fit to Equation (6) for the combined data 
set, the self-similarity of the two quite distinct pressure 
supported systems is evident; once scaled, both observed 
velocity dispersion data samples are consistent with each 
other to within their respective confidence intervals. Ac¬ 
tual physical radial scales vary from = 16.7pc of the 
globular cluster shown to Ra = 4.3A:pc of the elliptical 
galaxy. We see pressure supported systems as showing 
a simple empirical kinematic profile across many orders 
of magnitude in scale and distinct classes of systems. 

Despite the overlap evident in Figure 3, typically we 
find velocity dispersion profiles for elliptical galaxies to 
be rather flat, while those of Galactic globular clusters 
tend to be steeper. This is illustrated in Figure 4, which 
is analogous to Figure 3, but shows more typical ex¬ 
amples for the two astrophysical systems analyzed, the 
Galactic globular cluster NGC 7078 and the elliptical 
galaxy NGC 6515. Although the same functional form 
is evident, the degree of central concentration clearly 
varies as mentioned above. 

Regarding Galactic globular clusters, it has been sug¬ 
gested that the flattening of the velocity dispersion pro¬ 
files could be due to standard Newtonian physics, in 


particular tidal heating from the overall Milky Way po¬ 
tential (Kiipper 2010). Fortunately, 10 of the globular 
clusters included in our study form part of a wider sam¬ 
ple for which not only distances but also well measured 
proper motions are available. Indeed, Allen et al. (2006) 
and Allen et al. (2008) calculated orbital properties for 
a large sample of 54 Galactic globular clusters within 
a fully consistent Newtonian galactic model including 
baryonic disk, bulge and stellar halo, as well as a dom¬ 
inant dark matter halo. In those studies, the authors 
calculate detailed tidal radii for all globular clusters, not 
under any ’effective mass’ approximation, but through 
fully calculating the derivative of the total Galactic grav¬ 
itational force, including evaluating the gradients in ac¬ 
celeration across the extent of the clusters, along the 
entire orbital trajectories. 

In Hernandez et al. (2013a) one of us used fits to 
Equation (6) to obtain R„ for Galactic globular clus¬ 
ters including the ones in our present study, with the 
exception of NGC 2419, to compare with the Newto¬ 
nian tidal radii of Allen et al. (2006) and Allen et al. 
(2008). For the globular clusters in our present study 
(excluding NGC 2419) the ratio of the Newtonian tidal 
radii - at peri-Galacticon - to R^ ranges from between 
3 to 14.7, with a mean value of 7.1. Thus, it appears 
extremely unlikely than Newtonian tides are responsible 
for the flattening observed. The case of NGC 2419 is also 
clear, proper motions have only recently become avail¬ 
able due to the extreme distance at which it finds itself, 
Massari et al. (2016) provided the first proper motion 
studies for this cluster, and orbital integration within 
a Newtonian Galactic potential to show its galactocen- 
tric radius oscillates between 53 and 98 kpc, making it 
highly insensitive to Galactic tides. In going to MON- 
Dian schemes, as the gravitational force falls with R~^ 
rather than the gradient of it is much reduced, and 
Galactic globular clusters become more robust to Galac¬ 
tic tides than under Newtonian dynamics (e.g. Hernan¬ 
dez & Jimenez 2012), justifying their treatment here as 
isolated systems. 

For the case of the elliptical galaxies in the sample, we 
use two different criteria to determine interaction and 
isolation. The first comes directly from the morphologi¬ 
cal classification of the GALIFA 3rd data release, which 
is essentially a visual examination of the optical images 
of each object made by five different GALIFA members, 
where signs of mergers or interactions were sought. All 
our sample galaxies except one, NGG 4841A, appear as 
isolated under this criterion. In a second test, the / 
parameter discussed in Varela et al. (2004) is calcu¬ 
lated for all GALIFA galaxies (Walcher et al. 2014), 
which is an estimate of the logarithmic ratio between 
inner and tidal forces acting upon each galaxy, given 
in terms of the projected distance between target and 
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neighbor galaxies, the size of the target galaxy and the 
ratio of their masses as estimated from total luminosi¬ 
ties assuming a constant mass-to-light ratio. Galaxies 
with / values lower than —4.5 are considered completely 
isolated, while galaxies with / < — 2 are considered as 
not affected by tidal forces of their neighboring galax¬ 
ies. Again, all the galaxies from our sample except for 
NGC 4841A show / values lower than —2, indicating no 
tidal effects are evident upon them. In the case of NGC 
4841A, which shows clear visual signs of interaction and 
has a value of / = —0.01, a distinct perturbation in its 
velocity dispersion profile observations can be observed 
at around 12-14 kpc. The fit shown includes this slightly 
perturbed region, indicating the robustness of the fit to 
details; excluding all data beyond 12 kpc yields fit pa¬ 
rameters well within the reported confidence intervals of 
those resulting from using the full data sample. 

Other Newtonian explanations for the flattening in the 
globular cluster velocity dispersion profiles have been 
put forward in the literature, e.g. Kennedy (2014) 
showed that internal N-body relaxation processes can 
yield results consistent with observations. However, un¬ 
der any Newtonian explanation, it would then be down 
to a curious coincidence that the amplitude of CToo should 
scale with the total mass of the clusters in question 
(as inferred from detailed population synthesis models 
tuned to the particular ages and metallicities of each 
from McLaughlin & van der Marel 2005) in consistency 
with MONDian ’’Tully-Fisher” expectations. Equation 

(3), as shown in Hernandez et al. (2013a). 

We end this section with Figure 5, which contains the 
principal result of our study. The Figure shows Ra and 
(Too values for the fits to both of the pressure supported 
systems treated, in a logarithmic plane. The systems 
shown cover a range of 1.4— 165.6fcm/s in central veloc¬ 
ity dispersion and 1.8pc— 9.7kpc in the scale radius R^- 
The two distinct classes of objects are evident, with the 
small Galactic globular clusters appearing at the low (Too 
extreme as circles and triangles giving results for ellip¬ 
tical galaxies. The solid line is not a fit to the data but 
gives the Rm = 3CToo/ao prediction of MONDian gravity 
of Equation (4). 

In going from volumetric velocity dispersion profiles 
to the observable projected velocity dispersion profiles, 
various projection effects intervene, in ways which de¬ 
pend on the details of the degree of central concentra¬ 
tion of the volumetric velocity dispersion profile, and 
the real space density profiles of the tracers being ana¬ 
lyzed (Hernandez & Jimenez 2012; Jimenez et al. 2013a; 
Tortora 2014). Such details of the integration implicit 
in obtaining the line of sight projected profiles we treat 
necessarily introduce changing systematics across both 
classes of systems, from the centrally peaked kinematics 
of globular clusters to the flatter ones of elliptical galax¬ 


ies. To the above we must add the systematics inherent 
to any comparison across such distinct astrophysical sys¬ 
tems over such a range of distances and observational 
techniques as what we are attempting here. We thus 
find it encouraging of the physical interpretation pre¬ 
sented that the overall trend evident in Figure 5 can be 
clearly well represented by the predictions of Equation 

(4) , under the natural identification of Ra = Rm-, the 
transition to the flat velocity dispersion regime being 
given by the length scale of MONDian gravity. 

We also performed a weighted linear fit to the data in 
Figure 5 to equation; 

Such a linear fit to the full data set of Figure 5, shown 
by the dashed line, yields A = 0.73 ± 0.13 and n = 
1.89 ±0.32, in clear accordance with the MONDian pre¬ 
dictions of Equation (5) for A = 0.81 and n = 2. The 
corresponding ht to only the elliptical galaxies is largely 
unconstrained, but if we consider only the globular clus¬ 
ters we obtain Ago = 1.54±2.35and ncc = 1.45±1.12. 
This is not in conflict with the predictions of Equation 
5, although the overall trend is clearly driven by the po¬ 
sitions of the two clouds of points coming from the two 
distinct scale limited classes considered. 

Under Newtonian gravity, the flattening observed in 
elliptical velocity dispersion profiles, where tides are not 
relevant, would be interpreted as evidence for a dark 
matter halo, and the consistency of Ra oc cr^ of Figure 

(5) , as another curious coincidence down to the details 
of biasing and stellar feedback. 

Within MOND as such, for systems under the influ¬ 
ence of an external gravitational field, a regime change is 
expected according to whether the internal acceleration 
of the system in question is smaller or larger than the 
external one. In fact, if the external acceleration dom¬ 
inates and is larger than oq, the system being treated 
will behave in an essentially Newtonian manner, even if 
internal accelerations are well below qq (Famaey & Mc- 
Gaugh 2012). While not relevant for elliptical galaxies, 
the collection of Galactic globular clusters treated here 
span the above MOND transition, some of the closer 
ones would be expected to be affected by this so called 
external held effect, while the more distant, e.g. NGG 
2419, would be exempt from it, Sanders (2012). From 
Figure 5 it is clear that all clusters are consistent with 
the hrst order MONDian expectations of Equation (4), 
regardless of their galactocentric distances. Our results 
hence appear to support generic MONDian schemes 
where the external held effect might not appear, or re¬ 
sult in milder distortions than what appear in MOND 
as such, e.g. Verinde (2016) or Barrientos & Mendoza 
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R. Durazo, X. Hernandez, B. Cervantes Sodi & S. F. Sanchez 



log R [pc] 


Figure 5. (Too vs. Rrj values for both of the systems treated, 
red circles for Galactic globular clusters and green triangles 
for low rotation elliptical galaxies, in a logarithmic plane. 
The solid line is not a fit to the data, it rather shows the 
MONDian expectations of Equation (4) for the predicted 
scaling of Rm ~ 3a^/ao, Rm/ pc = 0.81{aoo/kms~^)^. 
A linear regression fit (dashed line) gives i?o- = (0.73 ± 
0.13 )(t^®®^°'®^, clearly compatible with the theoretical ex¬ 
pectation under the identification of Ra = Rm- 

(2016). 

It appears that a unique gravitational physics applies, 
with a clear transition radius at a Rm scale where the 
decline of the inner Newtonian R < Rm region gives 
way to the “Tully-Fisher” cr(i?) oc (GMoq)^/^ of the 
outer MONDian region, regardless of the fact that we 
start with self-gravitating systems of stars at pc scales, 
to the much larger ellipticals at kpc scales. 

Increasing our sample size and extending the sam¬ 
ple with other type of pressure supported systems, e.g. 
dSphs, to validate our conclusions is evidently a desir¬ 
able development, which however must be undertaken 
with care to keep the number of objects at each class 
considered relatively balanced, to make sure that the 


systematics inherent to any particular class are not bi¬ 
asing conclusions. Similarly, filling the gap present at 
(Too values of between 10 and 40 km/s would be highly 
desirable. 

5. CONCLUSIONS 

A universal projected velocity dispersion profile of 
(t{R) = + 0-00 has been shown to accurately 

model two distinct classes of astronomical systems; glob¬ 
ular clusters and low rotation elliptical galaxies. 

MONDian gravity generically predicts Rm = Scr^/oo, 
Rm/pc = 0.81((Too/fcms“^)^. Data across 7 orders of 
magnitude in mass give Ra/pc = A((Too/fcms“^)", with 
A = 0.73 ± 0.13 and n = 1.89 ± 0.32. 

Thus, in spite of significant observational uncertain¬ 
ties and varying projection effects across apparently very 
different classes of astronomical systems, the identifica¬ 
tion of the local physical Rm radius with the empirical 
projected kinematic R^ parameter yields results consis¬ 
tent with generic modified gravity expectations, where 
gravitational anomalies in the a < uq regime stem from 
changes in the physics and not a hypothetical dark mat¬ 
ter component. 
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